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In  this  report  we  present  recent  results  from  the  analysis 
of  lunar  laser  ranging  (LLR)  and  very  long  baseline  interfero¬ 
metry  (VLBI)  data. 

I .  Determinations  of  UTO  and  Variation  of  Latitude  from  Lunar 
Laser  Ranging  Data 


There  now  exists  over  ten  years  of  lunar  laser  ranging  data 
from  observations  made  by  the  McDonald  Observatory  near  Fort 
Davis,  Texas.  These  data  are  currently  being  analyzed  to  esti¬ 
mate  variations  in  the  rotation  of  the  earth.  The  primary  earth 
rotation  parameters  are  UTl  and  the  X  and  Y  components  of  pole 
position.  However  the  earth  rotation  parameters  to  which  LLR 
observations  are  fundamentally  sensitive  are  UTO  and  variation  of 
latitude.  In  this  respect,  the  information  content  of  LLR  obser¬ 
vations  from  a  single  station  is  the  same  as  that  of  observations 
from  a  photographic  zenith  tube.  In  order  to  determine  the  pri¬ 
mary  earth  rotation  parameters,  it  is  necessary  to  combine  data 
from  two  or  more  stations.  McDonald  Observatory  is  currently  the 
only  station  regularly  providing  scientifically  useful  data. 
Therefore,  if  we  ignore  other  sources  of  data,  we  must  content 
ourselves  with  determining  UTO  and  the  variation  of  latitude 
specifically  for  McDonald  Observatory.  On  the  other  hand,  if  an 
accurate  set  of  pole  position  values  is  available  from  a  separate 
source,  it  is  possible  to  use  them  in  conjunction  with  values  of 
UTO  at  McDonald  Observatory  to  determine  UTl.  Similarly,  if  a 


set  of  UTl  values  is  available,  it  is  possible  to  derive  the  X 
and  Y  components  of  pole  position  from  the  McDonald  Observatory 
UTO  and  variation  of  latitude  values. 

our  initial  approach  in  determining  earth  rotation  values 
from  the  LLR  data  was  to  determine  only  UTO  and  by  using  pole 
position  data  from  the  traditional  astronomical  techniques, 
determine  UTl.  We  realized  that  the  quality  of  the  pole  position 
data  used  would,  in  part,  determine  the  quality  of  the  UTl  es¬ 
timates.  In  this  regard  it  is  worth  noting  that  recently  signi¬ 
ficant  differences  have  been  found  between  the  values  of  pole 
position  and  UTl  obtained  using  traditional  techniques  as  re¬ 
ported  by  the  Bureau  International  de  I'Heure  (PIH)  and  those 
obtained  from  other  sources  (Stolz  e^  »  1976;  Harris  and 
Williams,  1977;  King  ^  1978;  Robertson  ^  Bl."  •  1980;  Bureau 

International  de  I'Heure,  1979].  The  discovery  of  these  differ¬ 
ences  has  resulted  in  the  establishment  of  a  new  BIH  system,  the 
"1979  BIH  System"  [Bureau  International  de  I'Heure,  1979].  The 
1979  BIH  System  was  adopted  in  an  attempt  to  "correct"  the  BIH 
results  to  bring  them  into  better  agreement  with  those  obtained 
from  the  other  sources.  Results  in  this  system  are  related  to 
results  in  the  1968  BIH  System  by  the  following  relationships: 

^1979  ~  ^1968  ""  ~  43932)] 

+  o"o07  sin[-~^~(MJD  -  43979)] 

M979  ^19r.8  ' 


UT1j979  ■  “'^^1968  "  0?0007  sin[^^|^(MJD  -  44048)] 

+  0f0OO7  sin[y~^{MJD  -  44019)] 

The  discovery  of  the  "errors"  in  the  BIH  results  shows  the  impor¬ 
tance  of  developing  and  refining  alternative  techniques  for 
determining  pole  position  and  UTl. 

Our  first  step  in  determining  UTO  was  to  use  the  MIT 
Planetary  Ephemeris  Program  (PEP)  to  estimate  tabular  values  of  a 
continuous  piece-wi se- 1  inear  model  of  UTO  with  respect  to  a  set 
of  a  priori  reference  values.  The  model  has  the  form 

^i+1  ' 

f(t)  =  g.  +  - ~(t  -  t.)  , 

i+1  i 

1  ^  5.  *-1  +  1'  i  =  » 

where  g^  =  f(tj)  are  the  unknown  tabular  points  which  are  esti¬ 
mated  simultaneously  with  all  other  parameters  affecting  the 
observations.  The  time  intervals  between  successive  tabular 
points  need  not  be  uniform;  they  were  chosen  at  approximately 
monthly  intervals  to  nearly  coincide  with  the  center  of  the 
usable  LLR  data  for  each  lunation.  The  number  of  tabular  points 
was  96. 

The  continuous  piece-wise-linear  model  provides  corrections 


to  the  reference  values  of  UTO.  These  corrections  contain  no 
information  on  variations  with  periods  less  than  about  one  month. 
In  a  subsequent  step,  we  estimated  the  high  frequency  variations 
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in  (JTO.  We  took  the  postfit  range  residuals  from  the  first  step 
and  analyzed  them  separately  for  one  observing  day  at  a  time  to 
obtain  an  incremental  correction  to  the  estimate  of  UTO  for  each 
day.  We  estimated  two  parameters,  a  range  bias  and  a  correction 
to  UTO  for  each  day  on  which  there  were  two  or  more  observations 
of  a  single  reflector  spanning  a  period  of  1.5  hours  or  more. 

Raw  estimates  of  corrections  to  UTO  were  obtained  for  636  dates 
between  October  1970  and  October  1979.  The  raw  daily  estimates 
were  smoothed  by  convolving  the  weighted  raw  values  with  a 
Gaussian-shaped  smoothing  window: 


-  x.)^/a^ 


*1 "  ‘ 
j=l 


where  the  weights,  w^,  are  given  by 


w 


j 


n 


i  =  l 


,2-2 

Xj)  /a 


-1 


with  the  (T^  being  the  standard  deviations  of  the  raw  estimates. 
The  smoothing  parameter,  _a,  determines  the  degree  of  smoothing 
obtained.  it  is  proportional  to  the  f ul 1-width-at-hal f-max imum 
(FWHM)  of  the  Gaussian  window: 


FWFIM  ==  1.66  5  a. 


Filters  with  FWHM  of  about  ^.1,  8.3  and  16.7  days  have  been  used. 

We  interpolated  among  the  smoothed  values,  using  cubic  poly¬ 
nomials,  to  obtain  estimates  of  UTO  at  5-day  intervals  corres- 
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ponding  to  the  dates  o£  the  published  BIH  values.  We  then 
combined  our  UTO  values  with  BIH  values  of  pole  position  to 
create  UTl .  A  comparison  of  these  values  with  smoothed  BIH 
values  in  the  1979  system  has  revealed  significant  differences 
while  a  comparison  with  a  limited  number  of  preliminary  VLBI 
(D.  S.  Robertson,  personal  communication,  1980)  and  raw  BIH 
determinations  showed  similar  trends  (see  Figure  1). 

The  accuracy  of  our  UTl  results  derive  partly  from  the 
assumption  that  the  pole  position  values  we  had  used  were 
accurate.  We  decided  to  test  this  assumption. 

The  histograms  of  Figure  2  illustrate  the  effects  of  nominal 
errors  of  1  ms  in  UTO  and  0"015  in  variation  of  latitude  on  the 
set  of  LLR  observations  through  May  1978.  The  units  of  delta  tau 
are  nanoseconds.  It  had  been  felt  for  some  time  that  the  inabi¬ 
lity  to  accurately  model  pole  position  might  be  significantly 
degrading  model  fits  to  the  LLR  data.  This  feeling  was 
dramatically  confirmed  by  a  simple  test  based  on  the  histogram 
results.  Using  the  full  set  of  observations  and  a  nominal  model 
which  included  BIH  values  of  pole  position  and  a  continuous 
piece-wise-linear  correction  to  BIH  values  of  UTl,  we  obtained 
postfit  residuals  which  had  a  normalized  RMS^  of  1.79.  Using 

1.  The  normalized  RMS  is  defined  as 

(  ^  (e./oj) ^/n) 
i  =  l 

where  n  is  the  number  of  residuals,  e^  is  the  value  of  the  ith 
residual,  and  a.  is  the  a  priori  estimated  standard  deviation  of 
the  ith  observation. 


MOmriED  JULIAN  DATE 

Fig.  1.  DIFFERENCES  BETOEEN  1)  LUNAR  LASER  RANGING  (LLR) ,  2)  VERY  LONG 
BASELINE  INTERFEROMETRY  (VLBI)  AND  3)  RAW  BUREAU  INTERNATIONAL  DE  L ’ HEURE 
(BIH)  VALUES  OF  UTl  AND  SMOOTHED  BIH  VALUES  OF  UTl  IN  THE  1979  SYSTEM. 
MODIFIED  JULIAN  DATE  43020  IS  28  AUG.  1976. 
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THE  SENSITIVITY  OF  TWO-WAY  RANGE  OBSERVATIONS  TO  ERRORS 
TO  AND  THE  VARIATION  IN  LATITUDE  DUE  TO  POLAR  MOTION 
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only  observations  for  which  an  error  of  1  ms  in  UTO  would 
introduce  an  error  of  less  than  0.2  ns  in  two  way  range  with  the 
same  model  as  used  for  the  full  set  of  observat ions ,  we  obtained 
postfit  residuals  with  a  normalized  RMS  of  1.61.  331  normal 

points  were  used.  Using  only  observations  for  which  an  error  of 
0"015  in  variation  of  latitude  would  introduce  an  error  of  less 
than  0.3  ns  in  two  way  range  with  again  essentially  the  same 
model  that  we  used  for  the  full  set  of  observations,  we  obtained 
postfit  residuals  with  a  normalized  RMS  of  only  1.28.  313  normal 

points  were  used  in  this  case.  These  results  indicated  that 
while  we  were  adequately  modeling  the  effects  of  variations  in 
UTO  on  the  LLR  observations,  the  variations  in  latitude  were  not 
adequately  modeled.  We  concluded  that  since  none  of  the  readily 
available  sets  of  pole  position  data  (BIH,  IPMS  f International 
Polar  Motion  Service] ,  DMAHTC  [Defense  Mapping  Agency 
Hydroqraphic/Topographic  Center])  are  sufficiently  accurate,  we 
should  estimate  the  variation  in  latitude  at  McDonald  Observatory 
directly  from  the  LLR  observations  themselves.  To  do  so,  we 
modeled  the  variation  of  latitude  with  a  continuous,  piece-wise- 
linear  functit)n  of  atomic  time  of  the  same  form  as  that  used  to 
determine  the  long  period  trends  in  UTO.  Again,  the  tabular 
point.s  were  estimated  simultaneously  with  all  other  parameters. 

The  time  intervals  between  successive  tabular  points  were 
the  same  as  those  used  for  the  function  describing  UTO.  The 
number  of  tabular  points  was  96.  By  adding  these  parameters  to 
our  analysis,  we  reduced  the  RMS  of  the  postfit  residuals  from 
2  6. 8  cm  to  19.6  cm,  a  reduction  of  almost  26%'. 
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We  have  compared  our  LLR  estimates  of  variation  of  latitude 
with  those  derived  from  the  pole  position  values  reported  by 
workers  using  other  techniques.  We  have  determined  the  variation 
of  latitude  at  McDonald  Observatory  from  (i)  the  smoothed  5-day 
Circular  D  values  of  the  BIH,  (ii)  the  smoothed  0.05-year  values 
of  the  IPMS  as  reported  in  the  Monthly  Notes  of  the  IPMS,  (iii) 
the  raw  daily  values  from  the  Polar  Motion  Reports  of  the  DMAHTC , 
(iv)  the  raw  5-day  values  produced  by  the  Goddard  Space  Flight 
Center  (GSFC)  from  laser  ranging  observations  of  the  satellite 
LAGEOS  and  (v)  the  raw  5-day  values  produced  by  the  Institute  of 
Advanced  Study  in  Orbital  Mechanics  (lASOM)  of  the  University  of 
Texas,  also  from  laser  ranging  observations  of  LAGEOS. 

The  raw  DMAHTC,  GSFC  and  lASOM  pole  position  values  contain 
a  significant  amount  of  high  frequency  variation.  Just  how  much 
of  the  high  frequency  component  is  noise  and  how  much  is  true 
pole  position  is  not  yet  known.  In  an  attempt  to  study  this 
problem  we  have  smoothed  the  pole  position  data  sets  and  produced 
smoothed  values  of  variation  of  latitude  at  McDonald  Observatory. 

Two  smoothing  techniques  have  been  investigated:  convolution 
of  the  raw  data  with  a  Gaussian-shaped  window  and  Whittaker's 
smoothing  procedure  [Whittaker  and  Robinson,  1940]  as  implemented 
by  Vondr^k  [Vondr^k,  1969;  Vondrak,  1977].  For  the  most  part  the 
procedure  of  applying  Gaussian  smoothing  to  the  pole  position 
data  follows  that  used  for  the  LLR  UTO  values.  However  for  the 
smoothing  of  the  pole  position  values  an  improvement  in  the 
Gaussian  smoothing  technique  was  made.  The  Gaussian  smoothing 
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technique  suffers  from  an  inability  to  provide  reasonable  values 
within  data  gaps  if  those  gaps  are  much  larger  than  the  full- 
width-at-hal f-max imum  (FWHM)  of  the  smoothing  filter.  This  is 
particularly  true  for  data  containing  periodicities  with  large 
amplitudes.  To  improve  the  smoothing  and  interpolating  ability 
of  the  Gaussian  smoothing  technique,  we  first  removed  from  the 
data  best-fitting  sinusoids  with  periods  of  432.1  days 
(Chandler),  365.24  days  (annual)  and  275.0  days  (a  periodicity  to 
be  discussed  later) .  We  then  smoothed  the  resulting  residuals 
and  subsequently  added  back  the  sinusoidal  terms. 

The  raw  DMAHTC,  GSFC  and  lASOM  values  were  smoothed  using 
Gaussian  windows  with  FWHM  of  about  8.3,  16.7,  33.3,  50.0,  66.6 
and  83.3  days.  The  variation  of  latitude  at  McDonald  Observatory 
was  determined  from  each  set  of  smoothed  values  and  compared  to 
the  estimates  derived  from  the  LLR  observations.  The  RMS 
differences  about  the  mean  were  calculated  for  each  comparison 
and  are  shown  in  Table  1. 
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Table  1 

RMS  DIFFERENCES  ABOUT  THE  MEAN  BETWEEN  LLR  AND  VARIOUS  OTHER 


DETERMINATIONS 

OF  THE 

VARIATION 

OF  LATITUDE  AT 

MCDONALD 

OBSERVATORY 

FWHM  (days) 

DMAHTC^ 

GSFC^ 

lASOM^ 

BIH68^ 

BIH79^ 

I  PMS^ 

8.3 

8.26 

16.92 

28.26 

16.7 

6.85 

14.51 

22.16 

33.3 

5.90 

11.23 

13.99 

8.57^ 

6.36^ 

1  1 .66^ 

50.0 

5.46 

9.49 

11.39 

66.6 

5.26 

83.3 

5.18 

1.  Organizations  are  defined  in  text.  Units  for  each  entry  are 
mil 1 iarcseconds . 

2.  BIH  and  IPMS  data  used  were  smoothed  by  those  organiza¬ 
tions.  The  equivalent  FWHM  of  a  Gaussian  filter  is  roughly 
one  month. 

The  time  span  used  for  ail  the  comparisons  was  MJD  =  ^2960  (10 
July  1976)  to  MJD  =  43809  (28  October  1978).  In  this  time  span 
each  data  set  contains  169  5-day  values.  The  overall  results  of 
this  comparison  had  been  anticipated.  Since  the  LLR  estimates 
are  in  effect  the  result  of  low-pass  filtering  with  a  cut-off  of 
roughly  one  month,  it  is  not  surprising  to  see  the  agreement 
between  the  LLR  estimates  of  variation  of  latitude  and  the  values 
derived  from  the  observations  of  other  techniques  improve  as 
these  values  are  more  heavily  smoothed.  Clearly,  the  best  agree¬ 
ment  between  the  LLR  values  and  those  of  the  other  techniques  is 
that  for  the  heavily  smoothed  satellite  Doppler  values.  However, 
if  the  raw  BIH  or  IPMS  data  were  to  be  as  heavily  smoothed,  they 
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miqht  also  provide  as  good  a  comparison.  We  intend  to  do  our  own 
smoothing  of  BIH  and  IPMS  data  to  confirm  this  prediction.  One 
cannot  judge  from  this  comparison  whether  any  of  the  high 
frequency  variations  seen  in  the  data  provided  by  the  other 
techniques  are  real. 

Another  comparison  of  pole  position  values  was  performed 
which  J_s  sensitive  to  high  frequency  variations.  Each  of  the 
data  sets  included  in  Table  1  (with  the  exception  of  BIH^S  and 
IPMS)  was  use(i  to  provide  a  priori  pole  position  values  in 
separate  analyses  of  tiie  I.LR  data  using  PEP.  Apart  from  the  pole 
position  values,  the  models  used  in  these  runs  were  identical. 

The  postfit  normalized  RMS  of  the  residuals  resulting  from  the 
analyses  are  shown  in  Table  2. 


Table  2 

POSTFIT  NORMALIZED  RMS  OF  LLR  RESIDUALS 
OBTAINED  USING  VARIOUS  SMOOTHED  POLE  POSITION  DATA  SETS 


FWHM  (days) 

DMAHTC 

GSFC 

lASOM 

BIH6B 

BIH79 

IPMS 

8.3 

1 .  38 

1 .83 

1  ,65 

ir,  .7 

1  .  30 

1.72 

1.56 

3  1.3 

1 . 28 

1.01 

1 . 80 

1 . 22 

SO  .0 

1  .22 

1.87 

1.49 

.f, 

1 . 20 

8  3.3 

1.10 

Again,  the  best  result  is  that  for  the  DMAHTr  data  smoothed  with 
a  filter  with  a  FWHM  of  83.3  days.  The  resiilt;-,  prest'nted  in 
Table  2  seem  to  indicate  that  a  lot  of  the  high  1  requenc'y  struc¬ 
ture  present  in  various  pole  positioti  data  set?;  is  noi-.e,  .'ith«‘r 
due  to  observat ional  error  or  reduction  error. 

The  Gaussian  smoothing  was  carried  out  after  the  best¬ 
fitting  Chandler,  annual  and  275-day  sinusoids  had  been  removed 
from  the  raw  data.  We  adopted  this  procedure  in  order  to 
improve  the  signal-to-noise  ratio  of  the  smoothed  data.  We 
observed  that  when  Gaussian  smoothing  was  carried  out  after 
sinusoid  removal  and  the  resulting  data  sets  used  to  provide  a 
priori  values,  there  was  a  drop  in  the  postfit  normalized  RMS  of 
the  LLR  residuals  over  that  obtained  when  the  pole  position  data 
were  filtered  without  sinusoid  removal.  As  expected,  the  drop 
was  particularly  significant  for  the  more  heavily  smoothed  data. 
For  example,  when  the  data  were  smoothed  with  a  filter  with  a 
FWHM  of  50.0  days  wi th  sinusoid  removal,  the  postfit  normalized 
RMS  of  the  LLR  residuals  was  10%  lower  than  the  normalized  RMS 
obtained  using  smoothed  data  wi thout  sinusoid  removal. 

All  of  the  DMAHTC ,  GSFC  and  lASOM  data  sets  used  in  the 
comparisons  presented  above  were  obtained  by  smoothing  the  raw 
data  (with  best-fitting  sinusoids  removed)  with  a  Gaussian 
filter.  The  other  smoothing  technique  which  we  have  investigated 
is  that  of  Whittaker.  This  work  is  being  carried  out  under 


another  contract  and  will  be  discussed  elsewhere. 
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In  the  course  of  removing  the  sinusoids  prior  to  smoothing 
the  DMAHTC  pole  position  values,  a  275-day  periodicity  in  the 
data  was  found.  The  approximate  period  of  this  term  was 
determined  from  a  visual  inspection  of  the  corresponding  peak  in 
the  frequency  spectrum  of  the  DMAHTC  data.  The  uncertainty  in 
the  estimated  period  has  yet  to  be  determined.  In  the  X- 
component  of  pole  position  this  275-day  term  has  an  amplitude  of 
13.1  mi  11 iarcseconds  (mas)  and  in  the  Y-component,  5.9  mas.  The 
phase  of  the  periodicity  in  the  Y-component  is  67°  ahead  of  that 
in  the  X-component.  This  term  has  been  observed  in  DMAHTC,  BIH  / 

and  IPMS  data  sets  by  other  workers  (Bowman  and  LeRoy,  1976; 

Graber,  1976;  Rudnick,  1956]  and  is  believed  to  be  real  although 
no  explanation  of  its  origin  has  yet  been  offered.  We  intend  to 
investigate  this  term  further. 

1 1 .  Performance  and  Analysis  of  Intercontinental  VLBI 
Exper iments 

We  have  conducted  a  number  of  intercontinental  VLBI  experi¬ 
ments  involving  the  Onsala  Space  Observatory  in  Sweden,  the 
Haystack  Observatory  in  Massachusetts,  and  several  other  United 
States  and  European  antennas.  The  results  from  experiments 
through  1978  were  analyzed  under  the  present  contract  and  have 
been  accepted  for  publication  by  the  Journal  of  Geophysical 
Research  (see  Appendix  A) . 

These  experiments  through  1978  were  carried  out  with  our 
Mark  I  VLBI  System.  Starting  in  November  1979,  our  interconti- 


nental  VLBI  experiments  have  utilized  our  new  Mark  iii  System 
which  is  more  than  tenfold  as  sensitive  as  the  Mark  I  System. 

Data  were  obtained  with  the  Mark  iii  System  in  experiments 
involving  Onsala  and  Haystack  in  November  1979,  July  1980,  and 
September  1980.  The  data  are  still  being  analyzed;  many  problems 
were  encountered,  as  described  in  part  in  our  Quarterly  Reports. 
However,  preliminary  results  are  now  available  from  the  November 
1979  experiment.  Table  3  shows  a  comparison  of  our  baseline 
estimates  from  these  data  with  our  earlier  determinations 
(Appendix  A).  The  agreement  is  reasonably  good,  the  consistency 
being  at  the  several-centimeter  level,  better  than  6  parts  in 

9 

10  .  Typical  postfit  residuals  for  the  group-delay  measurements 
are  shown  in  Figure  3  for  the  November  1979  VLBI  data;  the  RMS 
spread  is  only  about  30  ps,  equivalent  to  1  cm  in  distance. 
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GEODESY  BY  RADIO  INTERFEROMETRY; 

INTERCONTINENTAL  DISTANCE  DETERMINATIONS 
WITH  SUBDECIMETER  PRECISION 

T.  A.  Herring,*  B.  E.  Corey,*  C.  C.  Counselman  III,* 

I.  I.  Shapiro,*  ’’  B.  O.  ROnnSng,*  0.  E.  H.  Rydbeck,^ 

T.  A.  Clark,*  R.  J.  Coates, *C.  Ma , *  J.  W.  Ryan,* 

N.  R.  Vandenberg,*  H.  F.  Hi nteregger , * C .  A.  Knight,' 

A.  E.  E.  Rogers,"  A.  R.  Whitney, "D.  S.  Robertson,* 
and  B.  R.  Schupler* 

‘Department  of  Earth  and  Planetary  Sciences,  Massachusetts 
Institute  of  Technology,  Cambridge,  Massachusetts  02139. 

*Onsala  Space  Observatory,  S-439  00  Onsala,  Sweden. 

‘Goddard  Space  Flight  Center,  Greenbelt,  Maryland  20771. 
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Technology,  Cambridge,  Massachusetts  02139. 

Abstract .  Analysis  of  very-long-baseline  interferometer 
(VLBI)  observations  yielded  estimates  of  the  distances  between 
three  radio  telescopes  in  the  United  States  and  one  in  Sweden, 
with  formal  standard  errors  of  a  few  centimeters:  Westford, 
Massachusetts-Onsala ,  Sweden;  5, 599, 714. 66± 0. 03m;  Green  Bank, 
West  Vi rg ini a-Onsala ,  Sweden;  o , 31 9 , 31 7. 75± 0 .0 3m ;  and  Owens 
Valley,  Cal i fornia-Onsala ,  Sweden;  7 , 9 1 4 , 1 31 . 19± 0. 04m ,  where  the 
earth-fixed  reference  points  are  defined  in  each  case  with  re¬ 
spect  to  the  axes  of  the  telescopes.  The  actual  standard  errors 
are  difficult  to  estimate  reliably  but  are  probably  not  greater 
than  twice  the  formal  errors. 


1.  Introduction 


Very  long  baseline  interferometry  (VLBI)  has  been  used  for 
the  past  decade  to  determine  distances  between  radio  telescopes 
and  positions  of  radio  sources.  Improvement  in  the  precision  of 
baseline  length  determination  in  North  America  over  this  period 
has  been  substantial  --  from  a  precision  of  about  1  m 
[Hinteregger  et  al.,  1972)  to  a  precision  of  under  five  centi¬ 
meters  [Robertson  et  al.,  1979]  for  baselines  from  about  1000  km 
to  just  under  4000  km  in  length.  Near  the  middle  of  the  decade, 
distances  between  North  America  and  Europe  were  determined  with  a 
precision  of  about  half  a  meter  [Robertson,  1975;  see  also  Cannon 
et  al.,  1979].  In  this  paper  we  describe  more  recent  VLBI  deter¬ 
minations,  with  subdecimeter  precision,  of  these  intercontinental 
d  i  stances . 

2 .  Data  An  a  1  ys  i  s 

The  data  analyzed  consisted  of  interferometric  group  delays 
[Shapiro,  1976]  which  Vi/ere  obtained  from  three  sessions  of  obser¬ 
vations  involving  up  to  four  antennas,  as  described  in  Table  1. 
All  observations  were  made  with  the  Mark  I  VLBI  system  with  a 
center  radio  frequency  of  about  8  GHz.  A  multichannel  bandwidth 
synthesis  technique  [Whitney  et  al.,  1976]  was  used.  Each  re- 
Cfirded  channel  had  the  Mark  I  bandwidth  of  360  kHz,  whereas  the 
synthesized  tiandwidth  was  j-lOO  MHz  in  the  first  session  of  obser¬ 
vations  and  j'lOO  MHz  in  the  last  two  sessions.  A  hydrogen  maser 
frequency  standard  of  modern,  post-1970,  design  was  used  at  each 
telescope  for  each  session. 

The  Cartesian  coordinate  system  used  in  the  analysis  was 
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geocentric  and  earth-fixed  with  the  Z  axis  parallel  to  the  mean 
pole  of  rotation  of  1900-1906,  as  defined  by  the  International 
Latitude  Service  and  maintained  by  the  F5ureau  International  de 
I'Heure  (BIH)  .  The  X  axis  was  defined  to  be  pi'rpend  i  cul  ar  to  the 
Z  axis  and  in  the  direction  of  the  Greenwicli  meridian.  The  Y 
axis  completed  the  right-handed  triad.  Operationally,  the  origin 
of  this  system  was  defined  by  the  coordinates  for  the  intersec¬ 
tion  of  the  azimuth  and  elevation  axes  of  the  Haystack  radio 
telescope,  obtained  from  a  combination  of  spacecraf t-tracking  and 
VLBI  observations  made  at  various  sites.  The  orientation  was 
defined  by  the  BIH  values  for  pole  position  and  UTl  (1968  sys¬ 
tem),  with  the  addition  of  fortnightly  and  monthly  terms  of  small 
amplitude  [Woolard,  1959],  and  by  the  models  for  diurnal  polar 
motion  fMcClure,  19731  and  earth  tides  (based  on  ephemeris  posi¬ 
tions  of  the  moon  and  sun,  on  the  Love  numbers  h  0.61  and  !  = 
0.09,  and  on  the  assumption  of  no  dissipation).  The  small  ef¬ 
fects  of  antenna  deformations  (Carter  et  al.,  1980],  ocean  load¬ 
ing,  plate  tectonics,  and  other  errors  in  the  model  wore  ignored 
in  view  of  the  limited  precision  of  the  measurements  being 
analyzed. 

The  orientation  of  this  earth-fixed  coordinate  system  in 
inertial  space  (with  respect  to  the  mean  equinox  and  equator  of 
1950.0),  was  defined  by  the  pole  position,  by  UTl  and  by  the 
standard  formulas  for  sidereal  time,  nutation,  and  precession, 
with  certain  small  corrections:  Woolard's  [1953]  theory  as 
modified  by  Melchior  [1971]  was  used  for  nutation;  and  5,027?878 
per  tropical  century  at  1950.0  was  taken  for  the  precession 
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constant.  The  origin  of  right  ascension  was  defined  by  the  value 
12  hr  26  min  33.246  sec  for  the  radio  source  3C  273B  (elliptic 
aberration  removed) . 

Using  this  framework  and  a  theoretical  model  for  the  group 
delays  [Robertson,  1975;  Ma ,  1978],  we  estimated  by  weighted 
least  squares  the  baseline  vectors,  the  source  positions,  and  the 
other  relevant  parameters.  Included  in  the  last  category  are 
parameters  that  allow  for  the  estimation  of  changes  in  pole 
position  and  variations  in  UTl  when  data  from  two  or  more  ses¬ 
sions  of  observations  are  combined  [Shapiro  et  al.,  1974,  1976; 
Robertson  et  al.,  1979].  In  addition  to  these  purely  geometric 
parameters  there  are  parameters  in  the  model  that  represent  the 
effects  on  the  group  delays  of  the  propagation  medium  and  of  the 
behavior  of  the  clocks  at  the  various  radio  telescopes.  The 
representation  of  the  propagation  medium  includes  the  zenith 
t  ropospiier  i  c  delay  as  a  parameter.  Since  the  signal  delay 
through  the  troposphere  varies  with  time  and  with  location,  a 
number  of  such  parameters  are  utilized  v;ith  this  model.  The 
effects  of  the  ionosphere,  always  under  1  ns  in  the  zenith 
direction,  but  not  separately  modeled,  are  largely  absorbed  by 
these  parameters  because  of  the  similarity  in  the  signature  in 
the  group  delays  of  the  ionosphere  and  the  troposphere.  The 
representation  of  the  behavior  of  the  clock  at  one  site  relative 
to  tfiat  at  another  consists  of  a  polynomial  in  time  of  low  order 
with  the  coefficients  as  parameters,  the  first  two  denoting  the 
epoch  and  rate  offsets.  Since  the  clocks  drift  unpred ictably 
with  time,  a  number  of  such  polynomials  are  usually  employed. 
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The  choice  of  the  number  of  parameters  to  represent  the 
combined  effects  of  the  troposphere  and  the  ionosphere  and  the 
behavior  of  the  clocks  as  well  as  the  choice  of  the  time  interval 
of  applicability  of  each  parameter  are  to  a  certain  extent  sub¬ 
jective.  We  therefore  investigated  the  effects  on  the  results  of 
a  variety  of  such  choices.  As  extremes  for  the  representation  of 
the  troposphere  and  the  ionosphere,  we  considered  for  each  tele¬ 
scope  the  use  of  one  parameter  for  the  zenith  delay  for  each 
session  of  observations  and  one  such  parameter  for  each  12-hour 
period  of  observations.  For  the  clock  behavior  we  chose  the 
clock  at  Haystack  as  the  reference  and  considered  for  each  other 
telescope  for  each  session  the  use  of  one  parameter  each  for  the 
relative  offsets  in  epoch,  rate,  and  (in  several  instances)  rate 
of  change  of  rate.  The  intervals  of  applicability  of  such  'clock 
polynomials'  ranged  from  7.5  to  24  hours;  only  for  this  latter 
case  was  a  parameter  included  for  the  rate  of  change  of  rate.  The 
variations  in  the  estimates  of  baseline  lengths  resulting  from 
these  different  parameter i zat ions  were  up  to  1.8  and  1.2  times 
the  corresponding  standard  deviation  obtained  for  the  adopted 
parameterization  for  the  data  from  each  of  the  first  two  sessions 
and  from  the  third  session,  respectively.  In  view  of  this  sta¬ 
bility,  we  deem  it  unlikely  that  either  the  propagation  medium  or 
the  behavior  of  the  clocks  has  introduced  errors  in  baseline 
length  much  in  excess  of  the  formal  standard  deviations.  The 
adopted  choice  of  epochs  and  intervals  of  applicability  for  these 
parameters  (see  section  4)  was  based  primarily  on  examination  of 
the  postfit  residuals,  especially  those  from  observations  at  low 


-26- 


elevation  angles  for  the  parameters  representing  the  behavior  of 
the  propagation  medium. 

The  inverse  of  the  weight  for  each  group  delay  was  obtained 
from  the  sum  of  the  variance  found  from  a  signal-to-noise  analy¬ 
sis  [Whitney,  1974]  and  an  ad  hoc  term  included  to  account  for 
error  sources  that  were  not  a  function  of  signal  strength.  The 
inclusion  of  this  term  was  motivated  by  the  desire  to  {!)  allow 
for  the  effects  of  systematic  errors,  and,  relatedly,  (2)  weight 
the  observations  of  the  different  sources  more  evenly,  since  for 
some  strong  sources  the  signal-to-noise  analysis  yielded  uncer¬ 
tainties  that  were  well  below  the  suspected  contributions  of 
other  errors.  The  magnitude  of  this  ad  hoc  term  was  assumed 
constant  for  the  data  for  a  given  baseline  for  a  given  session  of 
observations.  Its  value  was  obtained  from  the  constraint  that 

I 

'  chi-square  per  degree  of  freedom  be  unity,  and  its  square  root 

ranged  from  0.10  to  0.25  ns.  Omitting  this  ad  hoc  term  from  the 
variances  changed  the  estimates  of  baseline  length  by  up  to  1.2 
times  the  standard  deviations  obtained  with  this  term  included. 

3.  Results 

Table  2  gives  our  estimates  of  the  distances  between  the 
radio  telescope  in  Sweden  and  the  three  in  the  United  States. 

The  reference  point  for  each  telescope  is  the  intersection  of  the 
azimuth  and  elevation  axes,  except  for  the  NRAO  telescope  where 
tfie  reference  is  the  point  on  the  polar  axis  closest  to  the 
(non i nter sect i nq )  equatorial  axis  [Hinteregger  et  al.,  1972]. 

Typical  postfit  residuals  from  the  simultaneous  analysis  of 
the  data  from  all  sessions  are  shown  in  Figure  1.  The  root  mean 
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square  of  these  residuals  ranged  from  0.13  ns  for  the  Haystack- 
NRAO  data  from  the  last  session  to  0.49  ns  for  the  Onsala-NRAO 
data  from  the  first  session,  the  ratio  being  approximately  as  ex¬ 
pected  in  view  of  (1)  the  threefold  smaller  synthesized  bandwidth 
used  in  the  first  session  of  observations,  (?)  the  generally 
larger  correlated  flux  densities  obtained  from  observations 
involving  the  shorter  baselines,  and  (3)  the  slightly  higher 
sensitivity  of  the  Haystack-NRAO  system  compared  to  the  Onsala- 
NRAO  system. 

The  repeatability  of  the  estimates  of  the  baseline  lengths 
shown  in  Table  2  indicates  that  the  formal  standard  errors  may 
not  be  much  smaller  than  the  actual  errors.  However,  one  fact 
detracts  from  the  importance  of  this  indication:  The  uncertain¬ 
ties  of  the  estimates  of  baseline  lengths  from  the  third  session 
of  observations  were  much  smaller  than  those  from  the  second 
which,  in  turn,  were  smaller  than  those  from  the  first.  The 
disparity  between  the  uncertainties  from  the  first  and  second 
sessions  was  due  primarily  to  the  contrast  in  synthesized  band- 
widths.  The  further  discrepancy  between  the  uncertainties  from 
the  first  two  sessions  and  those  from  the  third  was  due  mainly  to 
geometry;  the  schedule  of  observations  for  the  first  two  sessions 
was  controlled  by  requirements  of  other  experiments  involving 
only  the  U.S.  antennas  and  therefore  did  not  properly  provide  for 
Onsala's  participation. 

The  estimates  for  the  propagation  delay  through  the  atmos¬ 
phere  and  ionosphere  in  the  zenith  direction  ranged  from  6.7  ns 
(NRAO)  to  8.2  ns  (Onsala) ;  the  estimates  for  the  clock  parameters 
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were  also  each  in  approximate  accord  with  expectations,  ranging 
up  to  23  MS  for  epoch  offsets  and  from  >''3x10  to  «/'2xl0  '  for 

rate  offsets. 

The  estimates  for  the  source  positions,  the  distances  be¬ 
tween  the  radio  telescopes  in  the  United  States,  and  the  pole 
position  and  UTl  variations  will  be  given  elsewhere  and  discussed 
there  along  with  a  much  larger  collection  of  VLB!  observations 
that  involved  only  the  telescopes  in  the  United  States. 

4.  Conclusions 

Combination  of  data  from  three  sessions  of  VLBI  observations 
involving  sites  in  Sweden  and  the  United  States  yielded  results 
of  subdecimeter  precision  in  the  determination  of  the  inter¬ 
continental  distances.  Since  sufficiently  accurate  and  well- 
distributed  observations  to  determine  baselines  to  Sweden  with 
subdecimeter  precision  were  available  for  only  one  session,  the 
data  are  not  useful  for  the  establishment  of  subdecimeter  repeat¬ 
ability.  Neither  are  the  earlier  VLBI  determinations  [Robertson, 
1975]  useful  since  they  were  less  precise,  and,  in  addition, 
involved  a  different  telescope  at  Onsala  which  has  yet  to  be 
located  accurately  with  respect  to  the  one  we  utilized  (Table  1). 
Further,  to  assess  inherent  accuracy,  comparisons  would  best  be 
made  with  results  from  a  more  accurate  technique.  None  now 
exists.  The  most  accurate  alternate  technique  currently  avail¬ 
able  for  the  determination  of  these  particular  intercontinental 
baselines  is  based  on  observations  of  the  Doppler  shift  of  radio 
signals  transmitted  from  satellites,  plus  standard  ground  surveys 
to  tie  the  reference  points  for  the  Doppler  measurements  to  those 
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for  the  VLBI  determinations.  Such  Doppler  data  exist  (W.  E. 
Carter,  private  communication,  1979),  but  the  analysis  and  the 
ground  surveys  have  not  been  completed  and  so  no  comparisons  can 
yet  be  made. 

Periodic  repetitions  of  these  intercontinental  VLBI  measure¬ 
ments  are  also  planned  with  more  radio  telescopes  in  the  inter¬ 
ferometer  array  and  with  use  of  our  new,  more  powerful,  Mark  III 
VLBI  system.  This  system,  among  other  attributes  records  data 
simultaneously  from  two  well  separated  radio  frequency  bands, 
thus  allowing  ionospheric  effects  to  be  virtually  eliminated  as  a 
source  of  error.  The  purpose  of  these  further  VLBI  measurements 
will  not  only  be  to  check  on  the  repeatability  of  the  present 
baseline  determinations  but,  more  importantly,  to  improve  the 
precision  to  the  level  at  which  the  expected  [Minster  and  Jordan, 
1978]  '/’1. 7  cm/yr  spreading  rate  between  Europe  and  the  United 
States  could  be  detected  reliably  within  a  decade. 
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Adopted  solution 


-35  - 


Figure  Caption 

Figure  1,  Postfit  residuals  from  the  adopted  solution  (see  Table 

2)  for  the  Haystack-Onsala  group-delay  measuremenfs  from 
Session  3  (see  Table  1).  The  root-niean-square  of  these 
residuals  is  0.25  nsec.  The  vertical  bars  represent 
±  one  standard  deviation.  Each  letter  in  the  figure 
denotes  a  different  radio  source:  A  .  NRAO  150;  B  s 
3C  84;  C  =  4C  39.25;  D  ^  3C  273B;  E  ;  3C  345;  F  ; 

VRO  42.22.01;  G  '  2134+00;  H  -  3C  454.3;  and  I  H  OJ  287. 
One  other  source,  observed  on  some  baselines  in  Session 
3,  was  3C  279.  See  Clark  et  al .  (1976)  for  approximate 

values  of  the  coordinates  for  each  source. 


